A novel type of spiral wave with trapped ions 
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Pattern formation in ultra-cold quantum systems has recently received a great deal of attention. 
In this work, we investigate a two-dimensional model system accounting for the dynamics of trapped 
ions. We find a novel spiral wave which is rigidly rotating but with a peculiar core region in which 
adjacent ions oscillate in anti-phase. The formation of this novel spiral wave is ascribed to the 
novel excitability reported by Lee and Cross. The breakup of the novel spiral wave is probed and, 
especially, one extraordinary scenario of the disappearance of spiral wave caused by spontaneous 
expansion of the anti-phased core is unveiled. 
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Spiral waves are the most frequently encountered pat- 
tern formation in two-dimensional systems far away from 
equilibrium. They are usually thought to be responsible 
for the patterns in a wide range of systems (nonlinear 
optics [lj, magnetic films 0, new chemical systems 
subcellular Q biology, and complex plasma Q). The 
underlying dynamics supporting spiral waves are usu- 
ally oscillatory and excitable. Recently, spiral waves are 
also generated in media whose local dynamics is complex 
periodic or even chaotic @, 0]. Spiral waves can dis- 
play a number of distinct behaviors, some of which are 
quite complex. The simplest transition in spiral waves is 
Hopf bifurcation Q , which turns a rigidly rotating spiral 
wave into a quasi-periodic meandering one. Being an- 
other type of common transition, the breakup of spiral 
waves derives from two different ways: the core breakup 
and the far field breakup. In the core breakup, the spiral 
wave first develops into turbulence near the spiral core 
whose mechanism is due to the Doppler effects 0, [1(3] • In 
the far field breakup, the spiral wave first becomes un- 
stable far away from the spiral core and the instability 
underlying it is the absolute Eckhaus one [0, [H| ■ 

As one type of wave propagation, spiral waves are al- 
ways explored in systems where the spatial variable is 
continuous. However, spiral waves are also presented in 
systems of coupled oscillators where spatial variables are 
discrete. In the context of coupled oscillators, the states 
between adjacent oscillators are not required to be con- 
tinuous, which may give rise to some impressing phenom- 
ena on the dynamics of spiral waves. Kuramoto. el. al. 
studied a two-dimensional non-locally coupled oscillators 
and found the existence of spiral wave chimera 13|, [14 1 . 
In a spiral wave chimera, the oscillators in the core region 
of spiral wave are desynchronized while those around the 
periphery of the core are in synchronization. Martens et. 
al. analyzed the spiral wave chimera (l5| . Yang et. al. 
studied a two-dimensional locally coupled Rossler oscilla- 



tors |16| . They found a sandwiched spiral wave in which 
any two adjacent oscillators are in anti-phase and they 
attributed the presence of the sandwiched spiral wave to 
the shortwave instability of the homogeneous oscillation 
of the model [HI El]. 

Recently, pattern formation in ultra-cold quantum sys- 
tems which is of the nature of coupled system has recently 
received a great deal of attention. Lee and Cross consid- 
ered a chain of ions where dissipation is provided 
by laser heating and cooling and nonlinearity comes from 
trap anharmonicity and beam shaping. When the nonlin- 
earities and interaction are small perturbations relative 
to the harmonic motion of ions, they derived an ampli- 
tude equation for the ions 
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b is the coupling between adjacent trapped ions and c 
denotes how an ion's amplitude affects its harmonic fre- 
quency. The amplitude equation is similar to the com- 
plex Ginzburg-Landau equation(CGLE). Different from 
the CGLE which includes both reactive and dissipa- 
tive interactions, the above amplitude equation contains 
only reactive interaction since the adjacent ions interact 
through reactive Coulomb force. Lee and Cross consid- 
ered the pattern formation in the system described by 
Eq.(l); they found that long- wavelength waves are ex- 
pected when be > yet very short- wavelength waves for 
be < 0. In the case of no-flux boundary condition, the 
only allowed pattern formation for be > is homogeneous 
oscillation of all ions while very short-wavelength waves 
may display a more complicated structure for be < 0. 
The most surprising discovery by Lee and Cross is that, 
when the homogeneous oscillation of ions in the chain for 
6c > is perturbed by a localized pulse of anti-phase os- 
cillation, the perturbation would probably travel across 
the system for a long time before it dies off. This sort 
of excitability in an oscillatory medium is put down to 
the reactive interaction presented in the system. The 
new type of excitation of anti-phase perturbation to the 
homogeneous oscillation in Eq.(fT|) can also be observed 
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FIG. 1: (color online) 6=1 and c = 0.2. (a) The wave 
pattern of SWAPC for Re(Aij). (b) The snapshot of \A itj \. 
(c) The snapshot of (\Ai.j\)t- The middle four panels show 
the snapshots of Re(Aij) for different subsets and the bottom 
four show the snapshots of \ Aij \ for different subsets, (d) and 
(h) for the subset in which oscillators locate at positions with 
(2i,2j). (e) and (i) The subset with (2i, 2j + 1). (f) and (j) 
The subset with (2i + 1, 2j). (g) and (k) The subset with 
(2i + 1, 2j + 1). The color scale used by (a) [or (b)] are for all 
plots of Re(Aij) (or \Aij\) in this work. N = 50. 

when the system is generalized to a two-dimensional lat- 
tice. 

Consider that, for a spiral wave in a two-dimensional 
coupled system, its phase singularity (or the rotation 
center) may be off lattice and adjacent oscillators on 
the opposite sides of the phase singularity are always in 
anti-phase, which indicates that there exists a persistent 
source for anti-phase perturbation. Then it will be of par- 
ticular interest to delve how the new type of excitability 
reported by Lee and Cross influences the dynamics of 
spiral waves on a two-dimensional lattice. In this work, 
we will take this quest. 

We use a square lattice where ions are trapped to nodes 
on it. The system is described as 

dA- ■ 

= *6(-4A j|j + Aj-ij + A l+1J + Ay-i + A iJ+1 ) 

-(l + ic)|A 4J | 2 A 4J +Aij (1) 

with i,j = 1, N. Open boundary conditions as in Ref. 
flif are imposed upon the system. 

To numerically simulate Eq.([T]), we apply the fourth or- 
der Runge-Kutta method with dt = 0.001. A spiral wave 
occurs by initially setting the oscillators at the boundary 
with \Aij\ = 1 and the phases of Ajj running from 
to 2tt along the boundary. \Aij\ for other oscillators are 
set to be 0. We first let b =' 1, c = 0.2 and TV = 50. 
Stunningly, a novel type of spiral wave, a rigidly rotating 
spiral wave with a core in which adjacent oscillators are 
in anti-phase, shows up. 

For this novel type of spiral wave, two striking features 
can be revealed by the snapshots of the wave patterns of 



the real part of Ajj (Re(Ai.j)) in Fig. Ufa) and of the 
module of A(i,j) (\Aij\) in Fig. [Tfb). Firstly, the spiral 
wave possesses an odd core in which a pattern with the 
shortest wavelength appears [See Fig. [Ha)]; Secondly, 
there exists a transition annulus where the oscillation 
may be weaker (lower |j4(i,j)|). It is obvious that the 
transition annulus divides the wave pattern into the core 
region and the arm region [See Fig. [ljb)] and, in partic- 
ular, \Aij \ in the core region is as strong as that in the 
arm region. In Fig. QJc), we exhibit the pattern for the 
time average of \Aij\. Together with Fig. QJb), Fig. QJc) 
illustrates that the novel type of spiral wave is stable, 
confirmed by the unchanged location of the core. To be 
mentioned, the dimension of the core of the spiral wave 
may be defined as that of the transition annulus. Fur- 
thermore, we pick up two pairs of adjacent oscillators: 
one pair locates in the core region and the other in the 
arm region. The time evolutions of these four oscillators 
are monitored. As shown in Figs. EJa) and (b), all of 
them behave periodically. Nevertheless, it is interesting 
to find that the phase difference between adjacent oscilla- 
tors within the core region is around tt, which implies the 
anti-phase characteristic, and is also responsible for the 
pattern with the shortest wavelength in the core region. 
In contrast, there is only a minor phase shift between ad- 
jacent oscillators within the arm region, which is induced 
by the spiral wave propagation. In other words, oscilla- 
tors in the arm region perform an in-phased oscillation. 
In addition, Figs. [2ja) and (b) show that the anti-phase 
oscillating is greatly faster than that of the homogeneous 
one, which is also seen in Ref.0. Combining the above 
observations together, we draw that the novel spiral wave 
in Fig. [TJa) is a rigidly rotating spiral wave with an anti- 
phased spiral core (we denote it as SWAPC). 

In a CGLE, a rigidly rotating spiral wave always takes 
the following form: A(r,t) — F(r)exp{i[m6 + ip(r) — cot]} 
and F(r) = at the rotation center (or the phase sin- 
gularity) of the spiral wave 20]. In other systems such 
as reaction-diffusion systems, the similar formulation can 
be found provided that the spiral wave is rigidly rotat- 
ing. Whereas, the spiral wave presented in Fig. QJa) 
and (b) does not follow this formulation in the core re- 
gion. Particularly, the rotation center of this spiral wave 
is replaced by an anti-phased spiral core and the phase 
singularity of a normal spiral wave is lost. The statement 
is supported by the patterns of both Re(Ai.j) and \Aij\ 
for four subsets of the system: the subset consisting of 

all oscillators with location (2i, 2j)(i, j — 1, 2, , N/2), 

the subset with (2i + 1, 2j), the subset with (2i, 2j + 1), 
and the subset with (2i+l, 2j+l). As shown in Fig. QJd)- 
(g) , each subset displays a clear spiral wave pattern which 
shows that the anti-phased oscillation of Re(Aij) in the 
core region is modulated by the in-phased oscillation in 
the arm region. Nevertheless, \Aij\ in Fig. Qlh)-(k) f° r 
each subset substantiates that the phase singularity for 
a normal spiral wave is lost. Any subset itself possess 
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FIG. 2: (color online) The time evolutions of different oscil- 
lators locating at a line which crosses the core region of a 
SWAPC. 6=1 and c = 0.2. The subscripts (i, j) in these 
plots denote the positions of the oscillators, (a) A pair of 
oscillators in the arm region, (b) A pair of oscillators in the 
core region, (c)-(f) The successive oscillators in the transition 
annulus. The black curves are for Re(Ai,j) and the red ones 
for Im(Aij). (g) The evolutions of \Ai j\ for the oscillators 
in (c)-(f). ' 

no phase singularity characterized by \Aij\ = at the 
center of the spiral core in spite of the unique singularity 
a normal rigidly rotating spiral wave have. For each sub- 
set, two patches with low |A;j|, which situate oppositely 
to the center of the spiral core, turn up. It is no other 
than these patches in all subsets that contribute to the 
formation of the transition annulus. 

To further figure out the phase singularity for SWAPC, 
we regard how the in-phased oscillation in the arm region 
passes through the transition annulus to the anti-phased 
oscillation in the core region along a line extending trans- 
versely. We present the time evolutions of successive os- 
cillators on a line down the crossover area. The results 
given in Fig. EJc)-(f) show distinctly two periodic com- 
ponents in each oscillator. When the oscillator is close to 
the arm region, it is dominated by the in-phased oscilla- 
tion which is superimposed with a weak anti-phased oscil- 
lation. As the oscillator approaches the core region, the 
component of the anti-phased oscillation grows stronger 
and stronger and it becomes predominant. Recalling that 
the phase singularity of a normal spiral wave is actually 
manifested in the in-phased oscillation, it is the replace- 
ment of the in-phased oscillation by the anti-phased oscil- 
lation in the core region that causes the traditional phase 
singularity to be lost for SWAPC. To be mentioned, low 
\Aij\ (\Aij\ ~ 0) in the transition annulus is not re- 
lated to the phase singularity of a spiral wave. As seen 
from Fig. [2jg), low \Aij\ in the transition annulus just 
occurs when the anti-phased oscillation becomes compa- 
rable to the in-phased one. Especially, Fig. M,g) shows 
that \Aij \ ~ in the transition annulus emerges with an 
equivalent period as that of the anti-phased oscillation, 
which states that the existence of low \Aij\ in the tran- 
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FIG. 3: (color online) The spatial-temporal plots of Re(Aij) 
in (a) and \Aij \ in (b) for the development of SWAPC. b — 1 
and c = 0.2. 

sition annulus roots in the anti-phased oscillation and 
is foreign to the phase singularity. In short, the phase 
singularity in a normal spiral wave vanishes for SWAPC 
and, instead, an anti-phased core region plays the role of 
a rotor to support the spiral wave propagation. 

Emphatically, though the wave pattern in Fig. Ufa) 
looks like the spiral wave chimera, they are essentially 
different. Firstly, the phases of adjacent oscillators in the 
core region in a spiral wave chimera state are unrelated 
while the phase difference between adjacent oscillators 
in SWAPC is kept around it. Secondly, as discussed in 
Ref.[15|. the form A(r, t) = F(r)exp{i[m9 + ip(r) — u>t]} is 
recovered for a spiral wave chimera provided that A(r, t) 
is replaced by an order parameter R(r,t). That is, there 
exists a well defined phase singularity in the spiral wave 
chimera. However, it is quite different for SWAPC since 
the ordinary phase singularity has been substituted by 
an anti-phased core region where oscillation amplitude is 
as strong as that in the arm region. 

To get more insight into how the excitability of anti- 
phased perturbation to a homogeneous oscillation leads 
to the formation of SWAPC, we monitor the evolution of 
Re(Aij) and \Aij\ of a chain of oscillators crossing the 
spiral wave core from the very initial stage when the novel 
spiral wave begins to build. Both spatial-temporal plots 
of Re(Aij) and \Aij\ in Fig. [3] show a general scenario 
towards SWAPC. An ordinary spiral wave is generated 
originally; Then, near the phase singularity of the nor- 
mal spiral wave, an anti-phased perturbation comes into 
being; Ultimately, an anti-phased region with unchanged 
size appears and the novel spiral wave is established. The 
elaboration above explicitly demonstrates that the ex- 
citability of the anti-phased perturbation from the in- 
phased oscillation plays crucial roles on the formation of 
SWAPC. 

SWAPC could exist within a large parameter range 
when b and c alter. The variation of b or c affects on 
the dynamics of SWAPC distinctly. Augmented b always 
leads to the expansion of the anti-phased core and longer 
wavelength in the arm region. Oppositely, with fixed 
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FIG. 4: (color online) Different types of instability of SWAPC. 



the anti-phased perturbation brought about by the phase 
singularity, the normal spiral wave has to drift sponta- 
neously. Fig. Hie) shows a snapshot of a drifting spiral 
wave with normal phase singularity. The trajectory of 
the drifting spiral wave is presented in Fig. E{d) where 
the phase singularities at successive times are denoted 
on the plot. Clearly, the drifting spiral wave dies away in 
the end due to the collision between its phase singularity 
and the boundary. Then, we fix c while change b. For 
small b, one novel instability of SWAPC arises where the 
anti-phased core dilates spontaneously and persistently 
till overspreads the system and the arm region is even- 



The top panels show the snapshots of Re(Aij) and the bot- tually drive away. The snapshots of Re(Ai_j) and | ^4^ 



torn |i4i,j|. (a) and (b) are the far-away breakup of SWAPC 
owing to Eckhaus instability, b = 1 and c = —0.33. (c) and 
(d) are the evanesce of SWAPC resulting from spiral wave 
drifting. The green dots in (d) designates the trajectory of 
the successive phase singularities, b = 1 and c = 0.8. (e) 
and (f) are the novel near-core breakup of SWAPC due to 
the spontaneous expansion of the anti-phased core. 6 = 0.25 
and c = 0.2. (g) and (h) are the disappearance of SWAPC 
on account of the larger anti-phased core, b = 5 and c = 0.2. 
The size of the system N = 100. 



b, both the size of anti-phased core and wavelength in 
the arm region are independent of the variation of c. In 
addition, either an outwardly propagating SWAPC or an 
inwardly propagating SWAPC wave [2l| can be observed 
depending on b and c. 

To further investigate the dynamics of SWAPC, we 
investigate how it breaks up which is always an interest- 
ing issue in the field of pattern formation. We consider 
a large system with N = 100. The results are insen- 
sitive to the system size (We have verified this predic- 
tion for different system sizes). In a normal CGLE, the 
spiral wave always breaks up through Eckhaus instabil- 
ity. However, depending on how the parameters vary, 
SWAPC in the system Eq. [T] may display rich scenes. 
In this work, merely several simple situations are taken 
into account. We prescribe b = 1 at first. When c de- 
creases beyond a critical value, the Ekhaus instability 
steps in and the normal far-away breakup of spiral wave 
appears. A snapshot is shown in Fig. EJa) where a spi- 
ral wavelet with an anti-phased core is surrounded by a 
chaotic sea. The snapshot of the corresponding \A(i, j)\ 
presented in Fig. EJb) shows that there are many defects 
outside the spiral wavelet. As time goes on, the spiral 
wavelet will shrinks further and it will die off when the 
anti-phased core is swallowed by the defect sea. To be 
mentioned, the final state after the disappearance of the 
spiral wavelet with the anti-phased core is portrayed by 
the anti-phase pattern. On the other hand, though in- 
creasing c may also cause SWAPC to vanish, it does not 
induce spiral wave breakup. Actually, when c is above a 
critical value, SWAPC is replaced by a normal spiral wave 
in the beginning; furthermore, to suppress the growth of 



in Fig. He) and (f) exhibit a SWAPC with a bulk of 
anti-phased core. Different from aforementioned circum- 
stances, there is no instability of SWAPC is observed by 
increasing b. Nevertheless, there exist no SWAPC for suf- 
ficiently large b for any finite-size system since the ever 
swelling anti-phased core. 

In summary, we find a novel type of spiral wave when 
studying the dynamics of trapped ions on the model sys- 
tem Eq|T] a rigidly rotating spiral wave with an anti- 
phased core (SWAPC). The formation of SWAPC is 
mainly due to the excitability discovered by Lee and 
Cross. Despite the dynamics of SWAPC as well as its 
breakup scenarios has been thoroughly investigated in 
this work, there still remain some open problems such 
as how does the size of the anti-phased core change 
and what comprehensive description is contained in the 
breakup picture. Furthermore, whether SWAPC exists 
in other cold atom quantum systems is also a fascinating 
topic. 
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